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g : ABSTRACT 

We present high angular resolution observations with the Keck Interferometer, 
high dispersion spectroscopic observations with Keck/NIRSPEC, and near-IR 
photometric observations from PAIRITEL of a sample of 11 solar-type T Tauri 
stars in 9 systems. We use these observations to probe the circumstellar material 
within 1 AU of these young stars, measuring the circumstellar-to-stellar flux 
ratios and angular size scales of the 2.2 /im emission. Our sample spans a range of 
stellar luminosities and mass accretion rates, allowing investigation of potential 
correlations between inner disk properties and stellar or accretion properties. 
We suggest that the mechanism by which the dusty inner disk is truncated may 
depend on the accretion rate of the source; in objects with low accretion rates, the 
stellar magnetospheres may truncate the disks, while sublimation may truncate 
£C) [ dusty disks around sources with higher accretion rates. We have also included 

in our sample objects that are known to be highly variable (based on previous 
| photometric and spectroscopic observations), and for several sources, we obtained 

multiple epochs of spectroscopic and interferometric data, supplemented by near- 
q | IR photometric monitoring, to search for inner disk variability. While time- 

variable veilings and accretion rates are observed in some sources, no strong 
evidence for inner disk pulsation is found. 
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Subject headings: stars:pre-main sequence — stars: circumstellar matter — stars:individual(CI 
Tau, DK Tau, AA Tau, RW Aur, V1002 Sco, V1331 Cyg, DI Cep, BM And)- 
techniques:high angular resolution — techniques:interferometric 
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Introduction 



Protoplanetary disks are an integral part of the star and planet formation process, and 
studying such disks around early analogs of stars like our own sun can provide insight into 
how planetary systems form. T Tauri stars are < 2 M Q pre-main sequence objects, and their 
circumstellar environments are an ideal laboratory for studying the initial conditions of planet 
formation. A wealth of evidence, including direct imaging at millimeter and optical wave - 
lengths (e.g., iKoerner fc Sargentl Il995l ; iDutrey et al.l Il996l; iMcCaughrean &: O'Delll 119 96") . 



Adams et al. 



198 



Bertout et al- 



and modeling of spectral en ergy distributions (SEDs; e.g. 
19881 ; iBeckwith et al.lll990l ). has confirmed the long-espoused hypothesis that T Tauri stars 
are surrounded by substantial disks of dust and gas. Moreover, observed line profiles and 
UV continuum ex cesses indicate that T Tauri stars are accreting material from their circum- 



stella r disks (e.g., IWalkerlll972l ; iBasri fc Bertoutlll989l ; Edwards et al.lll994l ; iGullbring et al. 



19981) 



Study of the innermost disk regions may reveal the mechanism by which material is 
accreted through the disk onto the star. In the current paradigm, T Tauri disks are truncated 
by the stellar magnetosphere within the co-rotation radius, with material accreting alon g 
magnetic field lines onto high-latitude regions of the star (e.g.. lKonigllll99ll ; IShu et al.lll994j ). 
For typical T Tauri star masses, radii, magnetic field st rengths, and accre tion rates, predicted 



truncation radii range from ~ 0.02 — 0.2 AU (e.g., Ijohns-Krulll 120071 ). Previous studies 



have shown that the dusty c omponent of disks is generally truncate d further out than the 



magnetospheric radius (e.g.. Eisner et al.ll2005l ; lAkeson et al.ll2005bl ). Dust sublimation is 



caused by heating from the central star, plus possible additional heating due to viscous 
dissipation of accreting material and luminosity from the impact of the accretion flow onto 
the stellar photosphere. If the accretion luminosity varies with time, then the inner dust 
disk edge may pulsate as the sublimation radius moves in or out. Some e vidence for such 



variability may be present in near-infrared phot ometric monitoring data (e.g. JSkrutskie et al. 



19961 : ICarpenter et al.l l200ll ; Eiroa et al.l 120021 ) . Below, we discuss a search for inner disk 
pulsations using multiple epochs of spatially resolved observations at 2.2 fim. 

Knowledge of the spatial and temperature structure of protoplanetary disks at sub- 
AU radii is important for understanding the properties of dust and gas in regions where 
terrestrial planet formation occurs. Temperature is the crucial parameter in determining the 
radii where du st and water ice, key building blocks of terrestrial planets, can condense (e.g., 
Hayashilll98ll ). Spatial structure is important for u nderstanding how t he close-in extra-solar 
planets discovered by radial velocity surveys [e.g.. iMarcv et al.l 120051 1 either formed at, or 
migrated to, their observed orbital radii (e.g.. lLin et al.lll996l ). 



Infrared interferometric observations are necessary to spatially resolve disk regions 



- 3- 



within ~ 1 AU of the central stars. Such spatially resolved observations are the only way to 
probe directly the geometry and temperature of these inner disk regions. The Keck Interfer- 
ometer has recently enabled spatially resolved observations of protoplanetary disks around 



several approximate ly solar- type T Tauri stars (jColavita et al.ll2003l ; Eisner et al.ll2005l . 12006 



Akeson et al.ll2005al ); these observations constrained the dust disk inner radius and temper- 
ature, and showed evidence that the inner disk gas extends further in toward the central 
star. Moreover, when these data ar e combined with previous spatially resolved observations 



of more massive, luminou s objects (IMillan-Gabet et al 



2007; Leinert et al. 2004; van Boekel et al. 



2004 



199 



9. 



200 



Monnier et al. 



31; Eisner et al. 


2003. 


2004. 


20051. 


2006|; 


Akeson et al. 



2005bl: lLiu et alJl2007h. the inner disk properties as a function of luminosity can be inves- 



tigated (e.g.. IMillan-Gabet et al.ll2007l ). In general, these observations are consistent with 
expectations for disks with directly illuminated inner walls; however, there are notable ex- 
ceptions for solar and sub-solar mass T Tauri stars, motivating further observations. 

Here we present near-IR interferometric observations of a sample of low-mass (including 
solar and sub-solar mass) T Tauri stars. These observations spatially resolve the near-IR 
emission and allow measurements of the size of the emitting region. In addition, we use spec- 
troscopic and photometric measurements to constrain the relative stellar and circumstellar 
contributions to the near-IR emission, as well as the accretion luminosities of our targets. All 
of these data constrain the inner disk properties of our sample, and in cases where multiple 
epochs of data were obtained, enable study of inner disk variability. 



Observations 



2.1. Keck Interferometer Observations 



We used the Keck Interferometer (KI) to obtain spatially resolved measurements of the 
near-IR emission from a sample of 11 T Tauri stars (Tabled]). In addition, we incorporate 



i nto ou r dataset a measurement of RW Aur A obtained on 23 October 2002 by lAkeson et al. 



( I2005al ). Among the objects for which we have multiple epochs of data, DK Tau A, AA 
Tau, and R W Aur A are known to be photometrically variable at opti cal through near-IR 
wavelengths iHerbst et alill994l ; Iskrutskie et al.l fl996l : lEiroa et alil2002h . 



KI is a fringe-tracking long baseline near-IR Michelson interf erometer that combines 
the light from the two 10- m Keck apertures over an 85-m baseline (jColavita fc Wizinowich 
20031 ; IColavita et al.ll2003l ). For each target, we measured squared visibilities (V 2 ) at K- 
band (A = 2.2 //m, AA = 0.4 //m). The system visibility (i.e., the point source response 
of the interferometer) was measured using observations of unresolved calibrators (Table [1]), 
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weighted by the inter nal scatter in the c alibrator and the temporal and angular proximity 
to the target source ( Boden et al. 1998}) . Source and calibrator data were corrected for 



detection biases as described by lColavital (119991 ) and averaged into 5s blocks. The calibrated 
V 2 for a target source is the average of the 5s blocks in each integration, with uncertainties 
given by the quadrature addition of the internal scatter and the uncertainty in the system 
visibility. Systematic V 2 effects associated with enhanced performance of the AO system on 
one telescope relative to the other (which may d epend on the b rightness of the target) are 
calibrated by applying a "ratio correction" (e.g.. IColavitalll999l ). Typical uncertainties are 



-5%. 

Measured V 2 are shown in Figure [TJ The angular resolution of our observations is on 
the order of a milli-arcsecond (mas), which is sufficient to resolve the near-IR circumstellar 
emission, but not the central star. For expected stellar radii of 2-5 Rq and distances of 
140-700 pc for our sample, the stars subtend < 0.2 mas . Our measured V 2 thus include 
contributions from both the unresolved stellar component and the resolved circumstellar 
component. 



2.2. Keck/NIRSPEC Observations 

We obtained high- dispersion K-band spectra with NIRSPEC between November 2004 
and July 2006 for our sample (Table [p. We used the high-dispersion mode with the 3- 
pixel slit, which provides a resolving power of R ~ 24,000. The NIRSPEC-7 filter was 
used with an echelle position of 62.67 and a cross disperser position of 35.51. This provided 
seven spectral orders covering portions of the wavelength range between 1.99 and 2.39 /zm. 
Included in these orders are various stellar atomic lines, the CO (2-0), (3-1), (4-2), and (5-3) 
bandheads, and the Br 7 line. 



The spectra were calibrated and extracted using the "REDSPEC" package (e.g-. lMcLean et al. 



20031 ). Reduction included mapping of spatial distortions, spectral extraction, wavelength 
calibration, heliocentric radial velocity corrections, bias correction, flat fielding, and sky 
subtraction. We divided our target spectra by A0V or A1V stellar spectra (that were in- 
terpolated over broad Br7 absorption features) and multiplied by appropriate blackbody 
templates to calibrate the bandpass of the instrument and correct for telluric lines. For the 
two binary systems observed, which have separations of 2"3 (DK Tau AB) and l'.'A (RW Aur 
AB), the components are sufficiently well separated that their spectra could be extracted 
independently. 
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2.3. PAIRITEL JHK S Photometry 
We obse rved our sample wi th the 1.3-m Peters Automated InfraRed Imaging Telescope 



(PAIRITEL; iBloom et al.ll2006l ). and measured photometric fluxes at J, H, and K s bands. 
For each source, we obtained short exposures (52 ms) and relatively longer exposures (7.8 
s). Short exposures were read out as "clearing exposures" before the longer exposures (this 
is the nominal process of double-correlated reads in infrared imaging). The short exposures 
provide unsaturated images of our targets, while the longer exposures yield high signal- 
to-noise images of fainter stars within the observed fields that are used for photometric 
calibration. The camera provides 2" pixels over a 8'5 x 8'5 field of view. Images are obtained 
simultaneously in J, H, and K s as part of the camera design. 

Our data reduction consisted of background (sky+dark) subtraction (using dithered 
images to produce backgrounds) and flat fielding. Counts for target and calibrator objects 
were measured in apertures with 2-pixel radii, sky background was measured as the median 
value of pixels in an annulus spanning 10 to 20 pixels in radius, and the sky was subtracted 
from each pixel in the aperture. We calibrated the fluxes using measured 2MASS magnitudes 
for other stars within the field of view of each target image. 

Our photometric errors include scatter in measurements of source and calibrator counts, 
scatter in the magnitudes computed using different calibrators, and an assumed uncertainty 
in the 2MASS magnitudes of our calibrator stars of 0.017 mag. Measured fluxes and uncer- 
tainties for our sample are listed in Table [2j 

Where possible, we selected calibrator stars that did not appear to vary from epoch 
to epoch; since no photometric standard stars were observed, our photometry is measured 
relative to calibrator stars within the observed fields. We evaluated calibrator variability by 
comparison of each calibrator to the mean value of all calibrators within an epoch; as we 
discovered variable sources, we excluded them from the calculated mean value and iterated. 
However, for some fields where only a few stars were visible, no non-variable calibrator stars 
were found with this procedure (which may indicate unstable atmospheric or instrumental 
conditions); for these cases our photometric uncertainties include the apparent variability of 
the flux calibrators. 

The binaries in our sample (DK Tau AB and RW Aur AB) are not resolved with the 
2" pixels of the PAIRITEL camera. Thus the photometry for these objects is the composite 
of both components. The if -band brigh t ness ratios for DK Ta u AB and RW Aur A are 3.3 



and 4.3, respectively (ISimon et al.lll995l ; IWhite fc Gheal200ll ). and our photometry is thus 



dominated by the primary components. 
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3. Analysis 

In this section we use the data described in S|2] to constrain the properties of the inner 
disk regions around our sample objects. Our focus is on spatially resolving the circumstellar 
emission with our interferometry data. However, because we require measurements of the 
circumstellar-to-stellar flux ra tio in order to pro perly model the circumstellar component 



of our V measurements (e.g., lEisner et al.ll2005l ). we begin by deriving this ratio from our 



NIRSPEC data ( §3.ip . In §3.2[ we model our spatially resolved interferometric measurements 
of the circumstellar emission from our targets with a simple geometric model. In §3.3[ we 
derive accretion luminosities for our targets, to enable an investigation of how the inner disk 
geometries inferred from our interferometry data correlate with mass accretion rate. 



3.1. Circumstellar-to-Stellar Flux Ratios 

We use the NIRSPEC measurements described in §2.21 to derive 2 /im veilings for our 
sample. The veiling is defined as the ratio of the circumstellar excess flux to the stellar 
flux in the i^-band. Veilings are derived by comparing spectra of our target objects with 
spectra obtained (on the same night) for non-accreting, weak-lined T Tauri stars (WLTTs). 
These WLTTs provide excellent templates, since they have stellar effective temperatures 
and surface gravities similar to the classical T Tauri stars in our sample, but show neither 
near-IR excess emission nor accretion signatures. 

We select several atomic absorption lines for comparison between targets and templates. 
For a grid of WLTT templates, we perform a x 2 minimization between the lines seen in the 
target star and the same lines in the template stars after registration of the relative radial 
velocities, rotational broadening, and veiling. Radial velocity registration is accomplished via 
a cross-correlation of target and template star spectra. Rotational broadening is implemented 



by convolving the template spectrum with a Doppler line profile (lGraylll992l ). The template 



is veiled by adding a contin uum excess to the nor malized, un-veiled spectrum, and then 



re-normalizing the sum (e.g.. lBasri &: Batalhalll990l ): 



_ F u + r K 

r v, veiled 1 

l + r x 

For objects observed over multiple epochs, we first fitted r K and vsini simultaneously 
for each epoch, then fixed the v sinz as the weighted average of values from both epochs 
(which were consistent within uncertainties in all cases), and then repeated the fits with 
vsini fixed to this average value. Most of the fitted vsini values are smaller than the 
velocity resolution of our observations, ~ 12 km s~ x . This means that the vsini values of 
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these targets are compatible with those of the WLTT templates. Since the WLTTs used as 
templates in our veiling determinations may not be slow rotators, their spectra may already 
be rotationally broadened. Inferred values of v sin i for our targets thus represent only the 
amount by which the template spectra need to be further broadened, rather than the true 
v sin 2 of the source, and are only lower limits on the true vsini of our sample. 

For this fitting, we used Mg and Al lines between 2.10 and 2.12 /jm. We restricted 
the fits to 30 A regions around strong lines (Figure [2]) to avoid biasing the fits with noisy 
continuum regions of the spectra. As a check, we applied this fitting procedure to lines of Ti 
between 2.22 and 2.24 /im, in a different spectral order than the Mg and Al lines; the fitted 
veilings in both spectral orders agree within la uncertainties. 

When we fitted veiling to any of our WLTT templates, using the set of remaining 
WLTTs as templates for the fitting, we found veilings consistent with zero for all objects 
but one. For TYC 5882, the template with the earliest spectral type, the fitted veiling is 
0.2 ± 0.1; this is probably due to the mismatch between the spectral type of this and other 
templates (see Tabled]), and does not indicate that TYC 5882 is actually veiled. 

There may be some additional uncertainty in the inferred veilings for V1331 Cyg and 
DI Cep, since these sources have earlier spectral types (G5 and G8, respectively) than any of 
our templates (M2-K3; Table 1). To evaluate the effects of this source/template mismatch, 
we fitted our veiled, rotationally broadened te mplates to a non-accreting, ~ 30 Myr old 



G3/5 star (HIP 9141; IZuckerman fc Songi 12004 ) . We found that the best-fit veiling for this 



star is ~ 0.5, and we thus the veilings inferred for V1331 Cyg and DI Cep may be biased 
by this amount. We conservatively treat this bias as an uncertainty, added in quadrature to 
the statistical uncertainties in the fitted veilings; this additional uncertainty is included in 
Table EJ 

Our veiling measurements fall within ~ 1 month of our KI observations, providing 
estimates of the circumstellar-to-stellar flux ratios nearly contemporaneous with our inter- 
ferometric measurements. Moreover, this near-coincidence allows a comparison of veiling and 
the size scale of the circumstellar emission at multiple epochs (although there may be some 
variability in both quantities over the ~ 1 month timescales separating the spectroscopic and 
interferometric observations; see below). Unfortunately, no measurem ent of the if -band vei l- 



ing is available contemporaneous with the 2002 KI measurements of lAkeson et al.l (l2005al ). 
These authors estimated a circumstellar-to-stellar flux ratio of 0.6 by assuming that optical 
photometry traced the stellar photosphere, extrapolating the photosphere to if-band and 
then comparing the extrapolation to the measured if-band flux. Since the star may be veiled 
at optical wavelengths (due to emission from an accretion shock), the true photospheric flux 
may be lower than assumed, and their estimate of the if-band veiling is therefore a lower 
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limit. Nevertheless, we adopt the value of 0.6 here, and assume an uncertainty of 20%. 

Inferred veiling values are listed in Table [31 For A A Tau and RW Aur A, the veiling varies 
(at greater than the estimated la uncertainty level) from epoch to epoch. This indicates that 
the relative fluxes of the stellar and circumstellar components in these systems are changing 
with time. 



3.2. Inner Disk Sizes 



We now turn to analysis of our interferometric data. To determine the angular extent 
of the .fT-band emission, we analyze our KI data in the context of a circularly symmetric 
uniform ring mod el (e.g., Eisner et al.l 120031 ). This simple model has been shown (e.g., 
Eisner et al.l 120041 ) to provide a good description of more complex models of inner disks 
including puffed-up inner rims ( iDullemond et al.l 12001c llsella fc Nattall2005r). Further more, 
ring models have been used by previous investigators (e.g.. iMillan-Gabet et al.l 120071 ). and 
our use of this model will facilitate comparisons with previous results. 

We first rem ove the contr i bution to our measured V 2 from the unresolved central star. 
As described by lEisner et al.l ( 120051 ). decomposing the observed V 2 into stellar and cir- 



cumstellar components requires a measurement of the circumstellar-to-stellar flux ratio; the 
veilings determined in §3.11 provide this ratio for our sample. Contributions from the stellar 
and circumstellar components add linearly to the measured visibilities, and one can write, 



V 2 



disk 



y/^(l + F disk /^)-l 
-^disk/ F* 



(2) 



where F disk / is the circumstellar-to-stellar flux ratio. The V£ isk values for our sample, 
obtained using the veiling values measured closest in time to the measured V 2 , are plotted 
in Figure [31 



Uncertainties in our inferred circumstellar-to-stellar flux ratios lead to uncertainties in 
V^ sk in addition to those associated with V£ egs . This error is smaller when the circumstellar- 
to-stellar flux ratio is large. Given inferred veilings and associated uncertainties for our 
sample, the resulting additional uncert ainty is less than a fe w percent. For the first epoch 
of data from RW Aur A, adopted from lAkeson et al.l (j2005al ). the uncertainty is likely to be 
be larger ( §3.11) . 

When decomposing observed V 2 values into stellar and circumstellar parts, each epoch 
of KI observations was matched to the i^-band veiling measurement closest in time. Because 
our interferometric and spectroscopic observations are not exactly simultaneous (generally, 
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measurements are separated by a ~ 1 month), and objects may vary photometrically (at both 



optical and near-IR wav elengths) on this time-scale (e.g.. iHerbst et al.lll994l ; ISkrutskie et al. 



19961 ; lEiroa et al.ll2002l ). the inferred circumstellar-to-stellar flux ratios may not represent 



the true values at the time of the KI observations. Given this potential variability, we also 
decomposed our multi-epoch V^ eas into stellar and circumstellar components by adopting 
the average value of the veiling over all observed epochs. 

The circumstellar components of the visibilities, V 2 isk , were fitted with uniform ring 
models using a x 2 minimization to find the model providing the lowest residuals between 
predicted and measured V 2 . For DK Tau A,B, AA Tau, and RW Aur A, we fitted the ring 
model to each epoch of our KI observations. Best-fit radii for these face-on uniform ring 
models, using both the "nearest-in-time" and "average value" circumstellar-to-stellar flux 
ratio determinations, are listed in Table HI Uncertainties listed in the table are la statistical 



uncertainties in the fits, obtained by computing the surface where x 2 



A r 



+ 1. Both 



approaches to estimating the circumstellar-to-stellar flux ratios for multi-epoch observations 
yield sizes consistent within la uncertainties. In the analysis that follows, we will use the 
sizes derived with the nearest-in-time veiling measurements. 

Figure [3] shows that the V 2 for some objects (AA Tau, DK Tau A, RW Aur A, V1002 Sco, 
and AS 206) do not decrease with increasing projected baseline length; these data contrast 
with the predictions of circularly symmetric disk models. Photometrically determined stellar 
rotation periods combined with v sin % measure ments suggest highly inclined geometries for 
some of these sources (e.g.. lBouvier et al.lll995l ). consistent with the non-monotonic behavior 



of V 2 versus baseline length in Figure [31 However, because we are using a single-baseline 
interferometer (with a fixed separation and position angle), our data generally lack sufficient 
position angle coverage to place meaningful constraints on inclination (even accounting for 
Earth rotation, which helps to fill in the position angle coverage). 

If we fit our data (assuming that multipl e epochs of data t race a non-variable source) 
with an inclined uniform disk model (see e.g.. Eisner et al. 20031 ). we generally can not rule 
out an inclination of 0" (within 1-a nncertalntiS^ , we can not tell whether the non- 
monotonic behavior seen in Figure [3] is due to asymmetric source geometries or uncertainties 
in the measurements. For one source, RW Aur A, our data are sufficient to constrain the 
inclination. Assuming the source geometry has not varied over the two epochs of observation 
(see §4.2l for discussion of this assumption), the best-fit inclined disk model has an inclination 
of 77l}g degrees. 



1 For 3-parameter (inner disk radius, position angle, and inclination) fits, la uncertainties are computed 
from projections of the surface with x 2 = Xmin + 3.53. 
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3.3. Accretion Luminosities 



One of our aims in this analysis is to compare inferred inner disk properties with stellar 
and accretion properties. Accretion luminosity is a particularly interesting property to cor- 
relate with our inner disk measurements, since accretion can provide additional heating that 
modifies the inner disk structure. Here we use the Brackett gamma (Bry) emi ssion observed 
in our NIRSPEC data as a diagnostic of infall (IMuzerolle et al.l Il998l . 1200 ll ). and thereby 
estimate accretion luminosities. 

Br7 emission lines for our observed sources are shown in Figure HI We measure the Br7 
equivalent widths (EW) of our targets from these NIRSPEC observations by 1) determining 
the spectral baseline by fitting a straight line to a version of the spectrum that is interpolated 
over the Br7 feature; 2) removing this baseline from the (un-smoothed) target spectrum; and 
3) integrating under the Br7 line (over wavelengths from 2.157 to 2.175 /jm). We do not need 
to account for veiling here, since we are ultimately interested in the total line luminosity, 
which depends on the line EW as well as the continuum level (i.e., the correction to the EW 
due to veiling would be offset by de- veiling the continuum when computing the line flux). 

Uncertainties in these EW measurements arise from noise in the spectra and residual 
errors in calibration due to imperfect interpolation over the broad Br7 absorption features 
in the telluric standards. We assess these uncertainty sources by applying the same EW 
measurement procedure to non-accreting, weak-line T Tauri stars. We find that the la 
uncertainty level for EW estimates is 0.1-0.2 A for most sources. For one object, V1002 
Sco, the calibration did not work as well, and the error is 6 A. EWs and error bars for each 
source and epoch of observation are listed in Table [5j 

We use the Br7 EW measured from our NIRSPEC data and the A s -band magnitude 
determined from our photometry to estimate the Br7 line flux. For close binaries, where 
our near-IR photometry measures composite fluxes of the system, we es timate fluxes of in- 



dividual components using previously measured component flux ratios (ISimon et al.l 11995 



White fc Ghea l200ll ). which involves the implicit assumption that both component fluxes 
vary if the composite flux varie s. We de-redden our measured AVband magnitudes using 



A v values from the liter ature (IWhite fc Ghea 12001 



Eiroa et al. 



2002 



Steenman fc Thej ~ 1991 ). under the 



Walter et al. 1994 



as- 



Hamann &: Perssonlll992l ) and the reddening law of 
sumption that Ay is not time- variable. We convert these de-reddened magnitudes into units 
of flux density, F\. We then compute the flux in the Br7 line by multiplying F\ by the equiv- 
alent width: F Br7 = F\ x EW. The Br7 luminosity is then L Br7 = 4:iid 2 F Bn , where d is the 
dista nce. The Br7 luminosity is r elated to the accretion luminosity by an empirical relation 
from jMuzerolle et alJll99«L I2OO1I ): log(L acc /L ) = (1.26±0.19) log(L Br7 /^©) + (4.43±0.79). 
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Inferred values of L~B n and L acc are included in Table [5j Our es timated L acc values 



are compatible with previous measurements from the literature (e. g., 



Hartigan et al.lll995l : iGullbring et al.lll998l : ICalvet fc Gullbrindll998 



Valenti et a 



White k Ghez 



1993 



200lh . 



Previous measurements typically span an order of magnitude, and our inferred L acc are 
within an order of magnitude of all previously determined values. Moreover, the values for 
each star relative to others in our sample roughly match those of surveys in the literature; 
in surveys where several of our targets were observed, the lowest accretion rate objects are 
also the lowest accretors in that subset of our sample. 



4. Discussion 



4.1. Size-Luminosity Diagram 

If dusty disks extended inward to their central stars, their near-IR emission would be 
very compact, so much so that it would appear unresolved even with the mas-level angular 
resolution of near-IR interferometers. The fact that disks aroun d young stars are resolved 
with near-IR interferometric data (e.g.. iMillan-Gabet et al.lll999l ) indicates that these disks 
are truncated at stellocentric radii substantially larger than the stellar radiu s. Measured 
temperatur es at the disk truncation radii are typically ~ 1500-2000 K (e.g. . Eisner et al. 



2004. 120051 ). compatible with the sublimation temperature of silicate dust (e.g.. |Pollack et al. 



19941 ) . Determining how the inner disk truncation radius depends on source luminosity can 
constrain how t he inner edge is heated and, in turn, the radial an d vertical structure of the 



inner disk (e.g.. iMonnier fc Millan-Gabetl 120021 ; Eisner et al.ll2004l ). 



Figure [5] shows the inner ring radii determined for a sample of T Tauri and Herbig 
Ae/Be objects (including previous measurements from the literature) as a function of source 
luminosity. The source luminosity is the sum of the stellar luminosity and the accretion 
luminosity; in cases where no measurement of accretion luminosity is available (mostly high- 
mass stars), only the stell a r lum in osity is used. Inner ring si z es from the lite r ature are dra wn 
from IMillan-Gabet et al.l fepQlh : feisner et al.l J2004I . l2005h : kkeson et all feoool . l2005al jbh: 
Monnier et al.l J2005h : we used luminosities for these sources quoted i n these references. For 



the ob j ects in our current s a mple, we took stell a r luminosities from 



fll995h : IWhite fc Ghed fl200lf ); IWarter et all (119941 ) : lHamann & Persson 



Kenvon fc Hartmannl 



(1992). and accretion 



luminosities from Table [5j For objects where we measured the inner disk size over multiple 
epochs, the average size is plotted in Figure [5j 

We compare the size-luminosity diagram in Figure with the predictions of a simple 
physical model. Specifically, we assume that the inner disk is directly irradiated by the 
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central star, and that the inner edge therefore puffs up (iDullemond et al.ll200ll ). For this 
model, the location of the inner disk edge is, 



(1 + /) 



(3) 



Here, / is the ratio of the inner edge heig ht to its radius, and is e stimated to be between 0.1- 
0.2 for T Tauri and Herbig Ae/Be stars (IDullemond et al.ll200ll ). T; n is the temperature at 
the inner disk edge. Figure shows the inner disk radii predicted by this model for T in = 1500 
K and T in = 2000 K; these temperatures br acket the expected sublimation temperatures for 
silicate dust in a protoplanetary disk (e.g.. iPollack et al.lll994l ). 



Most objects with luminosities > 10 agree with the predictions of these mod- 



els; this fact has been no ted previously (IMonnier fc Millan-Gabetl 120021 ; lEisner et al.l 12004 



Millan-Gabet et al.l 120071 ). There are, however, several exceptional high luminosity sources 



for which geometrically thin disk models may be mor e suitable (e.g., lEisner et al.l 120041 ; 
Monnier fc Millan-Gabetl 120021 ; IVinkovic fc Jurkiq 120071 ). Moreover, Figure [5] shows that at 
the lowest stellar masses, this model may break down. The sizes of several low-mass sources 
are larger than predicted for these puffed-up inner disk models. 

We put forward several potential explanations for these deviations between puffed-up 
disk model predictions and data for low- mass T Tauri stars. First, the sublimation tem- 
perature may be reached further from the star. For example, viscous energy dissipation in 
accreting material may provide additional heating, pushing the sublimation radius outward. 
The sublimation temperature may also be lower, and hence located at larger stellocentric 
radii, in th ese low- mass objec ts because of systematically lower ambient gas densities (see 



Table 3 of IPollack et al.lll994l ). Small dust grains, which are hotter than larger grains at 



a given radius, may also lead to larger dust sublimation radii. Finally, these disks may be 
truncated by a mechanism other than dust sublimation, such as photoevaporative mass loss 
or magnetospheric truncation. In the following paragraphs, we argue that of the possible ex- 
planations considered, the most likely is a mechanism other than dust sublimation, namely, 
magnetospheric truncation. 

These various possible explanations make different predictions about the accretion rates 
of sources with exceptionally large inner disk radii. We therefore investigate how the discrep- 
ancy between measured sizes and predictions of irradiated disk models depends on a proxy of 
source accretion rat^E Figure [6] plots the difference between measured and predicted size as 



2 We use L acc rather than M, since L acc is closer to the observed quantities. Co nversion to M require s 



additional assumptions about stellar mass and radius, and the inner disk radius (e.g., Gullbring et al. 1998h . 
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a function of L acc /L*. We find that objects whose measured sizes are larger than predicted 
by models actually have lower (fractional) accretion luminosities than other sources. 



4-1.1. Accretional Heating 

If accretional heating causes larger inner disk truncation radii, then the low-mass stars in 
Figure [5] would have systematically higher mass accretion rates than other sources. However, 
Figure [6] shows that objects with large inner radii tend to have lower accretion rates. Thus, 
accretional heating seems an unlikely explanation for the large inner disk sizes. 



4-1.2. Lower Dust Sublimation Temperatures 

Dust sublimates more easily when ambient pressures are lower. Lower mass accretion 
rates lead to lower gas densities and pressures, and can therefore yield lower sublimation 
temperatures. The fact that larger disk truncation radii tend to occur around sources with 
lower accretion luminosities (Figure E]) is compatible with the hypothesis that these disks 
are truncated at lower dust sublimation temperatures. 

However, the dust sublimation temperature would need to be < 1000 K to explain 
many of the large measured disk sizes. A silicate dust sublimation temperature of 1000 K 
requires 8 orders of magnitude lower gas density than a sublimation temperature of 1500 K 



( iPollack et al.lll994l ). which corresponds to ~ 5 orders of magnitude lower mass accretion 
rates (assuming an a-disk to convert gas density into mass accretion rate). Such low accretion 
rates, and large dispersion in accretion rate, are not measured for our sample. It is therefore 
unlikely that lower sublimation temperature can explain the large truncation radii of disks 
around low-mass T Tauri stars. 



4-1.3. Small Dust Grains 



Smaller dust grains absorb and emit radiation less efficient ly than larger grains, an d 



this loss of efficiency is more severe at longer wavelengths (e.g., iBohren fc Huffman! Il983l ). 
As a result, small grains will achieve higher temperatures than larger grains at a given 
stellocentric radius. Thus, if a disk is composed of small grains, the dust sublimation radius 



(Monnier k Millan-Gabet 


2002; 


Isella et al. 


2006; 


2007 


). Using the results oi 


' Vinkovic et al. 


( 


2006 


)■ 
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our sample can be reproduced by puffed up inner disk models composed of grains sized 0.1- 
0.5 yum. While small-dust-grain models can reproduce these larger observed sizes, it seems 
unlikely that those objects that would require small dust grains to explain the observations 
are precisely those with the lowest mass accretion rates. 



4- 1-4- Photo-evaporation 



We now consider alternatives to dust sublimation, where the physical process respon- 
sible for truncating the disk might be more efficient at lower mass accretion rates. One 
such process is photo-evaporation, where the m ass-loss driven by s tellar UV radiation ca n 
overcome accretion at sufficiently low M (e.g., IClarke et al.l 1200 ll ; lAlexander et al.ll2006l ). 
The mass flow can "switch" from accretion to outflow when the accretion rate drops below 
~ 10~ 10 M© yr _1 . Since the objects in our sample have accretion rates > 10 -9 M yr _1 
(based on the accretion luminosities in Table EJ), photo-evaporation is unlikely to dominate 
accretion. Moreover, we do not see the signature of a discrete transition in Figure El at a cer- 
tain accretion luminosity. Thus, photo-evaporation is unlikely to be the physical mechanism 
behind the truncated disks observed for the low-mass T Tauri stars in our sample. 



4-1.5. Magneto spheric Truncation 



Large inner disk radii may also be due to magnetospheric truncation. T Tauri stars 
accrete material via magnetospheric accretion; the stellar magnetic field threads the disk and 
essentially creates a barrier to midplane accretion at the point of ba lance between magnetic 
pressure and inwa rd pressure f rom accretion (e.g., 
radius is given by (jKonigllll99ll ): 



Shu et al.l 119941 ) . The magnetospheric 



1?* 



2.27 



(S /l kG) 4 (i?*/Fl0 



1/7 



(4) 



,(M,/M )(M/1O- 7 M 

where -Bo is the dipole magnetic field strength. For many classical T Tauri stars, high disk 
accretion rates (> 10~ 7 M m yr -1 ) confine th is magnetospheric radius to well within the dust 
sublimation radius (e.g., Eisner et al.ll2005l ). and thus the observed inner edge of the dust 
disk occurs where temperatures reach the sublimation temperature (i.e., T m = 1500-2000 
K). 

However, lower accretion rates allow the magnetospheric radius to extend to larger 
stellocentric radii (-R mag oc M~ 2//7 ), while the lower associated accretion luminosity leads to 
a smaller stellocentric radius of dust sublimation. Low accretion rates can therefore cause 
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magnetospheric truncation outside of the dust sublimation radius (e.g.. Eisner et al.ll2006l ). 



If we take B = 1 kG, R* 
Tauri stars (e.g., 



2 R , and M* 



0.5 M Q as typical values for low- mass T 
Johns-Krulll 120071 ) . then mass accretion rates between 10~ 9 and 10 -8 M 



yr -1 (compatible with the lowest accretion luminosities in Table E]) yield magnetospheric 
radii between 0.07 and 0.15 AU, comparable to the measured inner disk sizes in Figure |5j 
We thus favor magnetospheric truncation as the explanation for the large inner disk sizes of 
low-mass T Tauri stars. 



4.2. Inner Disk Variability 



Many young stars, including those in our sample, are known to be photometrically 
variable at optical through near-infrared wayelengths, with variabi l ity timescales ranging 



from days to months (e.g.. iHerbst et al.l 11994 ; ISkrutskie et al.l Il996l ; ICarpenter et al.l 12001 



Eiroa et al.ll2002l ; iBouvier et al.ll2007l ) . While the underlying cause of the variability is usually 
attributed to star-spots, which modulate the emissio n on the stellar rotation period, in som e 



cases inner disk variability may also play a role (e.g., lEiroa et al.ll2002l ; IBouvier et al.l 120071 ). 



Time- varying increases in disk accretion rates, which might occur on inner disk dynam- 
ical timescales of days to weeks, could lead to increased disk heating and higher near-IR 
emission. Since a hotter disk could sublimate dust at larger stellocentric radii than a cooler 
disk, the inner edge of the dust disk may move in and out as the disk heating varies. Con- 
versely, if disks are magnetospherically truncated, then increases in accretion rate would lead 
to smaller disk sizes (Equation HJ). Variability in the stellar magnetic field strength can also 
cause variability in mass accretion rates and changes in the inner disk ra dius over timesc ales 
comparable to stellar dynamo oscillation period, typically ~ 10 yr (e.g.. lArmitagelll995l ). 



We obtained multiple epochs of spectroscopic and photometric data for DK Tau A and 
B, AA Tau, and RW Aur A and B, and multiple epochs of interferometric data for all of 
these except RW Aur B. As indicated in Tables [2HS1 several sources have variable infrared 
fluxes, veilings, and inferred accretion luminosities. Moreover, one object (RW Aur A) may 
have a variable inner disk size (Figure [3J Table H]). 

In the scenario where observed variability is due to changes in the disk heating (e.g., 
because of time- variable disk accretion rates), we would expect the following for higher disk 
temperatures: 1) i^-band flux will increase; 2) veiling will increase, since the circumstellar 
flux is increasing relative to the stellar flux; 3) accretion luminosity may increase (if this is 
the source of the disk heating); 4) the measured inner disk size should increase, since the 
dust sublimation radius would move out. 
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None of our targets clearly show such signatures. For DK Tau A, the f^-band veiling, 
near-IR fluxes, and inner disk size remain constant (within uncertainties) over the span of 
our multi-epoch observations. For AA Tau, the veiling appears to decrease from one epoch to 
the next while the Br7 flux also decreases. However, the inner disk size remains constant (or 
marginally increases; Table Hj)) and the -K" s -band flux appears to increase in the latter epoch, 
contrary to expectations for a dimmer disk. For RW Aur A, the Br7 luminosity appears to 
increase as the veiling decreases, which is not consistent with expectations for our scenario 
of inner disk variability above. DK Tau B and RW Aur B show no variability in any of our 
data. 

If the emission from our targets is circularly symmetric, then our modeling indicates 
potential variability (above the la level) in the inner disk size with time for RW Aur A. 
The inner disk radius appears to increase as the i^-band veiling increases, which is expected 
since a larger inner rim would produce more near-IR emission (for a fixed temperature). 
However, Figure [1] shows only marginal variability in the measured V 2 (in contrast to the 
V^ 2 isk plotted in Figure [3]) from one epoch to the next, and the apparent inner disk variability 
rests entirely on the inferred circumstellar-to-stellar flux ratios. The first epoch of KI data 
for RW Aur A does not have a reliable estimate of this ratio; the estimate of 0.6 from 



Akeson et al.l (l2005al ) is a lower limit ( §3.1|) . and a ratio as high as 2.5 would eliminate the 
apparent variability. Furthermore, a highly inclined disk can fit the data for RW Aur A well 
without any variability in the inner disk size, even if the inferred circumstellar-to-stellar flux 
ratios are taken at face value ( §3.21) . Thus we can not claim to see inner disk variability 
based on these data. 

The interpretation of our multi-epoch observations is complicated by the fact that dif- 
ferent data constituting each epoch were not obtained simultaneously. Data for a given 
epoch were obtained within a few week time-window. Because the dynamical timescales 
for inner disk variability are r dyn ~ a/ R 3 /GM < 10 days for the inner radii inferred for 
our sample, there may be substantial source variability between the different measurements 
within a single epoch. In future work, the observations constituting a single epoch should 
ideally be obtained within ~ 1 day of one another. In addition, a larger number of epochs 
would facilitate investigation of inner disk variability. 



5. Conclusions 

Using near-IR interferometric, spectroscopic, and photometric data, we measured the 
inner disk radii of 11 low-mass T Tauri stars. In addition, we measured the near-IR veilings 
and accretion luminosities of these objects. Our data substantiate previous claims that the 
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inner disk radii of T Tauri stars are generally consistent with the predictions of disk models 
with vertically extended inner rims at stellocentric radii where disk temperatures reach 
dust sublimation temperatures, between 1500 and 2000 K. However, there are a handful of 
sources, in particular objects with stellar luminosities < 1 L , for which the measured radii 
are substantially larger than predicted by such models. 

Discrepancies between models where disks are truncated at the dust sublimation radius 
and measured inner disk radii are larger for sources with low ratios of accretion to stellar 
luminosity. Thus, the physical mechanism by which the inner disk is truncated appears to 
depend on the relative importance of stellar irradiation and accretion. To explain this ob- 
servation, we suggest that at low accretion rates, as the pressure of accreting material drops, 
the point of balance between inward accretion pressure and outward stellar magnetic pres- 
sure moves to larger stellocentric radii. At sufficiently low accretion rates, magnetospheric 
radii can become larger than the dust sublimation radii, yielding the large observed sizes for 
low-accretion-rate objects in our sample. 

For several of our sample objects, we obtained multiple epochs of interferometric, spec- 
troscopic, and photometric data, and for each epoch the various datasets were obtained 
within a few weeks of one another. We used these multi-epoch observations to search for 
inner disk variability. While inferred veilings, accretion luminosities, and even inner disk 
truncation radii, appear to vary for some sources, none of the objects show a variability sig- 
nature consistent with changes in the inner disk. Future monitoring observations with more 
epochs, and with different data comprising each epoch obtained within a time-window short 
compared to inner disk dynamical timescales, will enable a more rigorous probe of inner disk 
variability. 

The near-IR interferometry data presented in this paper were obtained with the Keck 
Interferometer (KI), which is supported by NASA. We wish to thank the entire KI team 
for making these observations possible. KI and NIRSPEC observations were carried out at 
the W.M. Keck Observatory, which is operated as a scientific partnership among California 
Institute of Technology, the University of California, and NASA. The Observatory was made 
possible by the generous financial support of the W.M. Keck Foundation. The authors wish 
to recognize and acknowledge the cultural role and reverence that the summit of Mauna 
Kea has always had within the indigenous Hawaiian community. We are most fortunate 
to have the opportunity to conduct observations from this mountain. This work has used 
software from the Michelson Science Center at the California Institute of Technology. J.A.E. 
acknowledges support from a Miller Research Fellowship, and thanks G. Basri for useful 
discussions and input into various stages of this work. J.A.E. is also grateful to F. Ciesla for 
interesting discussion about dust sublimation temperatures, and to E. Lopez for his interest 
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Table 1. Target and Calibrator Properties and Observations 



Source 




5 


d 


Type 


KI 




NIRSPEC 


References 


Target Stars 


CI Tau 


04 33 52.01 


+22 50 30.1 


140 


K7 


24 Oct 2005 




07 Jan 2006 


1 


DK Tau A 


04 30 44.28 


+26 01 24.6 


140 


K9 


30 Oct 2004, 24 Oct 2005 




19 Nov 2004, 07 Jan 2006 


2 


DK Tau B 


04 30 44.42 


+26 01 23.9 


140 


Ml 


30 Oct 2004, 24 Oct 2005 




19 Nov 2004, 07 Jan 2006 


2 


AA Tau 


04 34 55.45 


+24 28 53.7 


140 


K7 


30 Oct 2004, 24 Oct 2005 




19 Nov 2004, 07 Jan 2006 


1 


RW Aur A 


05 07 49.568 


+30 24 05.161 


140 


KI 


30 Oct 2004* 




20 Nov 2004, 07 Jan 2006 


3 


RW Aur B 


05 07 49.465 


+30 24 04.82 


140 


K5 


30 Oct 2004 




20 Nov 2004, 07 Jan 2006 


3 


V1002 Sco 


16 12 40.505 


-18 59 28.14 


160 


K0 


21 Apr 2005 




01 Jun 2005 


4 


AS 206 


16 25 56.3 


-24 20 50.0 


160 


K6.5 


21 Apr 2005 




13 Jul 2006 


5 


V1331 Cyg 


21 01 09.21 


+50 21 44.8 


700 


G5 


24 Oct 2005 




13 Jul 2006 


6,7 


DI Cep 


22 56 11.534 


+ 58 40 01.795 


300 


G8 


24 Oct 2005 




13 Jul 2006 


5,8 


BM And 


23 37 38.47 


+48 24 12.0 


440 


K5 


24 Oct 2005 




07 Jan 2006 


9,10 


KI Calibrator Stars 


HD 28459 


04 30 38.408 


+32 27 27.965 


94 


B9.5V 


30 Oct 2004, 24 Oct 2005 






11 


HD 282230 


04 31 23.148 


+31 16 41.390 


77 


K0V 


30 Oct 2004, 24 Oct 2005 






11 


HD 30378 


04 48 22.743 


+29 46 22.811 


196 


B9.5V 


30 Oct 2004, 24 Oct 2005 






11 


HD 36724 


05 34 58.506 


+ 26 58 19.958 


86 


F7V 


30 Oct 2004, 24 Oct 2005 






11 


HDC142943 


15 57 57.772 


-17 05 23.974 


75 


F6V 


21 Apr 2005 






11 


HDC144925 


16 09 02.600 


-18 59 44.037 


150 


A0V 


21 Apr 2005 






11 


HD 214279 


22 35 51.994 


+56 04 14.067 


108 


A3V 


24 Oct 2005 






11 


HD 214240 


22 35 53.382 


+50 04 14.836 


500 


B3V 


24 Oct 2005 






11 


HD 219290 


23 14 14.385 


+ 50 37 04.408 


126 


A0V 


24 Oct 2005 






11 


HD 219891 


23 19 02.413 


+45 08 12.302 


126 


A5V 


24 Oct 2005 






11 


NIRSPEC Calibrator Stars 


HR 1251 


04 03 09.3801 


+05 59 21.498 


40 


A1V 






19 Nov 2004, 07 Jan 2006 


11 


HR 1252 


04 59 15.4096 


+37 53 24.881 


49 


A1V 






20 Nov 2004 


11 


HR 5511 


14 46 14.9241 


+01 53 34.388 


39 


A0V 






13 Jul 2006 


11 


HR 8518 


22 21 39.3754 


-01 23 14.393 


48 


A0V 






01 Jun 2005 


11 


NIRSPEC WLTT Template Stars 


AG Tri A 


02 27 29.2543 


+30 58 24.613 


42 


K6 




19, 


20 Nov 2004, 07 Jan 2006, 13 Jul 2006 


11 


AG Tri B 


02 27 29.2543 


+30 58 24.613 


42 


M0 




19, 


20 Nov 2004, 07 Jan 2006, 13 Jul 2006 


11 


TYC 5882 


04 02 16.5 


-15 21 30 


63 


K3/4 






20 Nov 2004 


12 


RX J0409.8 


04 09 51.108 


+24 46 21.20 


140 


Ml. 5 






19, 20 Nov 2004, 07 Jan 2006 


13 


RX J0432.8 


04 32 53.237 


+ 17 35 33.65 


140 


M2 






19, 20 Nov 2004, 07 Jan 2006 


13 


RX .10437. 4A 


04 37 26.89 


+ 18 51 25.9 


140 


K6 






19, 20 Nov 2004, 07 Jan 2006 


13 


RX .10437. 4B 


04 37 26.89 


+ 18 51 25.9 


140 


M0. 5 






19, 20 Nov 2004, 07 Jan 2006 


13 


RX J0438.2 


04 38 15.619 


+23 02 27.88 


140 


M2 






19, 20 Nov 2004, 07 Jan 2006 


13 


RX J1534.3 


15 34 23.1 


-33 00 07 


145 


M0 






13 Jul 2006 


14 


RX J1546.0 


15 46 05.5 


-29 20 40 


145 


M0 






01 Jun 2005, 13 Jul 2006 


14 


RX J1546.7 


15 46 47.0 


-32 10 06 


145 


M2 






01 Jun 2005 


14 


RX ,11557.8 


15 57 50.0 


-23 05 09 


145 


M0 






01 Jun 2005 


14 


RX J1558.5 


15 58 53.6 


-25 12 32 


145 


Ml 






01 Jun 2005 


14 


Note. — *-for RW Aur A, we 


also utilize KI mca 


surcmcnt 


s from 


23 Oct 2002 jAkeson et al.ll2005 


,). 


All Taurus sources were assumed to be at 


a distance of 



140 pc. 



References — (1 ) iKenvon fc HartmanrJ jl995f); (2) IMonin et aTl il998h ; (3) IWhite fc Ghezl )200lD; (4) IWalter et al.l tl994l) ; (5) iHerbig fc Belli il988l); (6 ) 
IChavarrialcl <1981l>; (7 ) lHamann fc Perssonl 11993): ( 81 IKholopov j (195911; (9) l Aveni fc Hunter! 11969 ); (10) iMora et al.l )200lT) ; (11) iPerrvman et al.l 1199711 : 
(12) IZuckerman fc Song ( 20041); (13) IWichmann et al.l J20001 ); f 14) iKohler et al.l J200ol ) 
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Table 2. PAIRITEL JHK S Photometry 



Source 


Date 


mj 


m H 








lau 


Z4UCLU0 


q t^n j- n no 
o.OU zt U.UZ 


7 OO -L n 

( .yy zt u 


no 
.UZ 


7 ci in 
( .01 zt U 


nQ 
.Uo 




ZUINOVUO 


o ok j- n no 
o.zo zt U.Uz 


(.00 zt U 


no 
.UZ 


7 Q/l -L n 

( .o4 zt U 


no 
.UZ 


\)\\ lau A 15 


18UCLU4 


7 c\a j_ n no 
/ .y4 zt U.Uz 


7 AR -L n 

( .40 zt U 


no 
.UZ 


7 nQ j- n 
( .Uo ± U 


nQ 
.Uo 




UD1NOVU4 


7 -U n no 
/ .yo zt u.uz 


7 A r; _i_ n 
* .40 zt U 


no 

.UZ 


a Cio _i_ n 
o.yo zt U 


no 
.UZ 




Z4UCtU0 


7 qq j- n nrc 
f .00 zt U.UO 


7 QQ -L n 

( .oo zt U 


n^ 
.Uo 


c\ no _L n 
o.yz zt u 


no 
.Uo 


A A Tom 

aa lau 


loUCtU4 


& &n -i- n no 

O.oU zt U.UZ 


o. Q7 -i- n 

O.O 1 ± U 


no 

.UZ 


7 Q/i _i_ n 
* .y4 zt u 


HQ 
.Uo 




U01NOVU4 


o ifi-Ln nQ 
y.iu zt u.uo 


o.Oo zt U 


nQ 
.Uo 


O.IZ zt U 


nQ 
.Uo 




z4UCtUD 


q co j- n nQ 
o.oZ zt U.Uo 


Q QQ -L n 

o.oo zt U 


nQ 
.Uo 


7 7/1 -L n 

i . / 4 zt U 


(\A 

.U4 




ZUINOVUO 


c fir; _i_ n nQ 
o.oo zt U.Uo 


o. on -i- n 
o.zu zt u 


nQ 
.Uo 


7 QO -1- n 
( .oZ zt U 


nQ 
.Uo 


T$XK1 Anv A 13 

xvW Aur A 15 


18UCLU4 


c ^q j- n no 
o.Oo zt U.UZ 


o in i n 

o.iy zt u 


no 
.UZ 


7 CtQ j_ n 
f .Do zt U 


no 
.UZ 




24Oct05 


8.51 ±0.02 


8.00 ±0 


.02 


7.54 ±0 


.02 


V1002 Sco 


21Apr05 


8.29 ±0.02 


7.87 ±0 


.02 


7.35 ±0 


.02 


AS 206* 


21Apr05 


9.0 ±0.3 


8.5 ±0 


.3 


7.9 ±0 


.3 


V1331 Cyg* 


12Sep04 


10.01 ±0.04 


9.64 ±0 


.04 


8.98 ±0 


.04 




13Oct04 


10.02 ±0.04 


9.65 ±0 


.04 


9.06 ±0 


.04 




23Oct04 


10.03 ±0.04 


9.57 ±0 


.04 


9.00 ±0 


.04 




17Nov05 


10.06 ±0.04 


9.63 ±0 


.04 


9.03 ±0 


.04 




24Oct05 


10.06 ±0.04 


9.67 ±0 


.04 


9.07 ±0 


.04 




25Oct05 


10.01 ±0.04 


9.64 ±0 


.04 


9.04 ±0 


.04 




26Oct05 


9.93 ±0.04 


9.54 ±0 


.04 


8.95 ±0 


.04 


DI Cep* 


13Oct04 


9.1 ±0.2 


8.6 ±0 


.2 


8.2 ±0 


.2 




20Oct04 


9.2 ±0.2 


8.8 ±0 


.2 


8.3 ±0 


.2 


BM And 


13Oct04 


9.85 ±0.03 


9.38 ±0 


.03 


8.89 ±0 


.03 




23Oct04 


10.12 ±0.03 


9.59 ±0 


.03 


9.14±0 


.03 




24Oct05 


9.61 ±0.02 


9.10 ±0 


.02 


8.61 ±0 


.02 




17Nov05 


9.84 ±0.03 


9.32 ±0 


.03 


8.85 ±0 


.03 



Note. — *-For AS 206, V1331 Cyg, and DI Cep, our photomet- 
ric calibrators appear to be variable at a low level, and the quoted 
uncertainties include this potential error. 
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Table 3. fT-band veiling measurements 



Source 


Date 




Rotational Broadening (km s 1 ) 


Template Star 


CI Tau 


07 Jan 2006 


1.16 ±0.13 


< 12 


RXJ 0437.4A 


DK Tau A 


19 Nov 2004 


1.58 ±0.16 


< 12 


RX J0437.4B 




07 Jan 2006 


1.32 ±0.14 


< 12 


RX J0437.4B 


DK Tau B 


19 Nov 2004 


0.78 ±0.10 


< 12 


RX J0437.4B 




07 Jan 2006 


0.80 ±0.15 


< 12 


RX J0437.4B 


AA Tau 


19 Nov 2004 


0.54 ±0.06 


< 12 


RX J0437.4B 




07 Jan 2006 


0.34 ±0.06 


< 12 


RX J0437.4B 


RW Aur A 


20 Nov 2004 


o 9fi +0.86 
°- zu -0.62 


35 ±10 


AG Tri A 




07 Jan 2006 


1.50 ±0.20 


35 ±10 


AG Tri A 


RW Aur B 


20 Nov 2004 


0.32 ±0.08 


11±4 


RX J0437.4B 




07 Jan 2006 


0.28 ±0.06 


11 ±4 


RX J0437.4B 


V1002 Sco 


01 Jun 2005 


0.36 ±0.08 


78 ±6 


RX J1557.8 


AS 206 


13 Jul 2006 


1.48 ±0.22 


8±7 


AG Tri A 


V1331 Cyg 


13 Jul 2006 


> 10 






DI Cep 


13 Jul 2006 


3.40 ±0.71 


9±5 


AG Tri A 


BM And 


07 Jan 2006 


O Qn+ - 44 
z - ou -1.50 


17 ±17 


TYC 5882 



Note. — The Tk and rotational broadening values were determined by minimizing the \ 2 
residuals between observed spectra and a grid of veiled, rotationally broadened template spectra 
(20J- The template stars are WLTTs showing no near-IR excess emission, and fitted veilings 
for these objects are consistent with zero. However, the templates may have non-negligible 
wsini values, and thus the fitted rotational broadenings for our targets are lower limits for the 
true vsini values. Quoted error bars are la statistical uncertainties in the best-fit parameters. 
Spectral types of template stars can be found in Table Q] 



-26- 



Table 4. Face-on Ring Radii 



Source 


Date 




, (AU) 


i? ring (fixed r K ) 


CI Tau 


24Oct05 


0.097 


± 


0.008 


0.097 


± 


0.008 


DK Tau A 


30Oct04 


0.103 


± 


0.005 


0.105 


± 


0.005 




24Oct05 


0.107 


± 


0.008 


0.105 


± 


0.008 


DK Tau B 


30Oct04 


0.137 


± 


0.018 


0.137 


± 


0.012 




24Oct05 


0.137 


± 


0.018 


0.138 


± 


0.018 


AA Tau 


30Oct04 


0.084 


± 


0.029 


0.090 


± 


0.031 




24Oct05 


0.133 


± 


0.019 


0.120 


± 


0.016 


RW Aur A 


23Oct02 


0.103 


± 


0.005 


0.111 


± 


0.007 




30Oct04 


0.181 


± 


0.016 


0.131 


± 


0.007 


RW Aur B 


24Oct05 


0.161 


± 


0.034 


0.161 


± 


0.034 


V1002 Sco 


21Apr05 


0.118 


± 


0.027 


0.118 


± 


0.027 


AS 206 


2lApr05 


0.112 


± 


0.007 


0.112 


± 


0.007 


V1331 Cyg 


24Oct05 


0.315 


± 


0.021 


0.315 


± 


0.021 


DI Cep 


24Oct05 


0.165 


± 


0.039 


0.165 


± 


0.039 


BM And 


24Oct05 


0.249 


± 


0.028 


0.249 


± 


0.028 



Note. — i?ring is computed from fits of uniform ring mod- 
els to KI data, and the distances assumed in Table Q] Quoted 
error bars are la statistical uncertainties in the % 2 minimza- 
tion between models and measured V^ 2 isk . Column 3 lists 
the ring radii determined when V^ 2 isk are computed using the 
veiling measurement closest in time to the KI measurement. 
Column 4 lists the ring radii for V^ 2 isk estimated using the 
average veiling measured over all epochs. 
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Table 5. Br7 Emission Properties 



Source 


Date 


EW (A) 


L Bn (10- 4 Lq) 


L acc (L©) 


L* (Lq) 


CI Tau 


07 Jan 2006 


7.4 ±0.1 


2.3 


0.70 


0.8 


DK Tau A 


19 Nov 2004 


7.6 ±0.1 


2.6 


0.80 


1.7 




07 Jan 2006 


6.9 ±0.1 


2.5 


0.76 


1.7 


DK Tau B 


19 Nov 2004 


2.7 ±0.1 


0.3 


0.05 


0.5 




07 Jan 2006 


1.8 ±0.1 


0.2 


0.03 


0.5 


AA Tau 


19 Nov 2004 


4.8 ±0.1 


0.7 


0.16 


0.6 




07 Jan 2006 


0.2 ±0.1 


< 0.1 


< 0.01 


0.6 


RW Aur A 


20 Nov 2004 


9.3 ±0.1 


1.5 


0.42 


1.7 




07 Jan 2006 


16.5 ±0.1 


3.1 


1.02 


1.7 


RW Aur B 


20 Nov 2004 


6.7 ±0.1 


0.3 


0.05 


0.4 




07 Jan 2006 


5.1 ±0.1 


0.3 


0.04 


0.4 


V1002 Sco 


01 Jun 2005 


14 ±6 


5.4 


2.07 


3.8 


AS 206 


13 Jul 2006 


9.7 ±0.2 


2.4 


0.74 


1.6 


V1331 Cyg 


13 Jul 2006 


14.4 ±0.2 


24.5 


13.78 


21.0 


DI Cep 


13 Jul 2006 


8.4 ±0.2 


4.6 


1.69 


5.2 


BM And 


07 Jan 2006 


1.4 ±0.1 


1.0 


0.23 


5.5 



Note. — Given the uncertainties in the photometry and EWs, we estimate uncer- 
tainties between 5% and 20% for Lr^. The conver sion from L-Qi-y to L acc includes 
additional uncertainties ( Muzerolle et al. 19981 2001 ). and L acc is thus uncertain by 
a factor of a few. We also list stellar lum i nosities, drawn from the literature 
J White fc Ghezl l200ll: Irlartigan et~aH Il995l : IWalter et all Il994l ; lHamann fc Persson 
1992b . 
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Fig. 1. — Measured V 2 from KI as a function of uv radius, r uv (r uv is equal to the projected 
baseline length divided by the wavelength). Given the angular resolution of our observations, 
these measurements include contributions from the unresolved central star and resolved 
circumstellar emission. Different epochs are indicated with different symbols: unfilled circles 
repre sent the second epoch (where available). The first epoch of KI data for RW Aur A are 
from lAkeson et al.l (j2005al ). 
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Fig. 2. — Portions of NIRSPEC spectra for R W Aur A (solid cur ves), with the wavelengths of 



photospheric Mg and Al absorption lines (e.g.. |Prato et al.ll2003l ) indicated. The spectrum of 



a non-accreting WLTT (AG Tri A), which has been rotationally broadened and veiled (and 
vertically offset from the target spectrum), is also plotted (dotted line). Spectral regions 
used when matching rotationally broadened, veiled template spectra to target spectra are 
indicated by straight horizonal lines. In 2004 Nov (left), photospheric absorption lines from 
Mg and Al were significantly more veiled than in 2006 Jan (right). The inferred veilings 
are Tk = 3.26 and 1.50, respectively. We note the presence of emission from the H2 line at 
2.12182 /im in both epochs; however analysis of this feature is beyond the scope of this work. 
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Fig. 3. — Circumstellar component of the measured V 2 from KI (V 2 isk as computed from 
Equation [2]). These squared visibilities represent the disk component only. Plotted error bars 
include uncertainties in measured V 2 and in circumstellar-to-stellar flux ratios. Different 
epochs are indicated with different symbols: unfilled circles repres ent the second epoc h 
(where available). The first epoch of KI data for RW Aur A are from lAkeson et al.l (j2005al ). 
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Fig. 4. — Br7 lines observed for our sample. For objects observed over multiple epochs, the 
line observed during the first epoch is plotted with a solid line, and the line observed during 
the second epoch is represented with a dotted line. 
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Fig. 5. — Inner ring radii determined from near-IR interferometry, as a function of source 
luminosity (stellar+accreion luminosities, where available). Inner ring sizes from the lit- 
erature (filled circles) and for the sources in our sample (open squares) are plotted. W e 



indic ate the inner ring radii expected for a puffed-up inner disk model (iDullemond et al. 



200ll ) with T in = 1500 K (dotted line) and T in = 2000 K (dashed line). While most objects 
with luminosities > 10 L agree with the predictions of these models (with the excpetion of 
several high-mass stars whose disks appear consistent with geometrically thin disk models; 
Eisner et al.l 12004 ). at lower luminosities, there are numerous sources with inner ring radii 
larger than predicted values. 
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Fig. 6. — Discrepancy between the measured size and the model prediction (in the range of 
inner disk temepratures between 1500 and 2000 K) is plotted against the ratio of accretion 
to stellar luminosity. The discrepancy is the difference between the measured size and the 
region bounded by the two curves in Figure El As in Figure El disrepancies for sources 
from the literature are indicated by filled circles and, for the sources in our sample, by 
open squares. We have restricted this plot to T Tauri stars (L* < 5 L ), since most higher 
luminosity objects do not have measured L acc . While some low accretion luminosity sources 
do agree with the model predictions, a trend is nevertheless evident whereby objects with 
smaller L acc /L* ratios are more discrepant from model predictions. The objects exhibiting 
the largest discrepan cies are DK Tau B and RW Aur B from our sample, GM Aur from 
Akeson et all J2005ah . and AS 207 A from feisner et all J2005h . 



